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1. Introduction 

Of the large variety of topological defects that are known, domain walls 
are the simplest to understand in many ways. The underlying topology 
that gives rise to domain walls is discrete and can easily be visualized. In 
some models the equations of motion can be solved analytically leading to 
explicit expressions. Since the topology is discrete, formation via the Kibble 
mechanism can be discussed in terms of known results in percolation theory. 
Hence it is instructive to embark on a study of topological defects starting 
with the simplest case of domain walls. The purpose of these lectures is to 
get you started. 

As often happens, the simplest of considerations evolve into more com¬ 
plex phenomena. And domain walls are no exception. Only in the simplest 
of models are domain walls simple. As we go to somewhat more elaborate 
and realistic models, the domain walls in them get more complicated and 
realistic! Some of the developments that I will describe have direct connec¬ 
tion with the investigation of domain walls in He-3 [1]. In these lectures 
I will develop the subject enough so that you can see the complexities 
emerging. However, there will not be enough time to cover all the recent 
advances and those woods will be left for you to explore on your own. An 
important part of the forest that I will not get into is the cosmology of 
domain walls. This includes discussion of their dynamics and cosmologi¬ 
cal evolution. There simply isn’t enough time for it. I am hoping that the 
lectures by Professor Arodz will rectify some of these omissions. 

I will begin by describing the simplest of domain walls. This is the “kink” 
in a Z 2 model. Then I will consider the more complicated field theoretic 
model based on SU (5) x Z 2 symmetry. (The discussion is easily carried over 
to the case of SU(N) x Z 2 for odd N.) The discrete symmetry responsible 
for domain walls is still Z 2 like in the simple model. Here, however, we will 
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find a discrete spectrum of domain walls all having the same topology but 
different masses and other properties. The formation of domain walls in this 
model via the Kibble mechanism and the importance of these considerations 
in connection with cosmology will be discussed. 

2. The kink 

Consider the 0 4 Lagrangian in 1+1 dimensions labeled by (t,z) 

L =\{ d »< t >) 2 - V 2 ) 2 ( 1 ) 

where (f>(t, z ) is a real scalar field - also called the order parameter. The 
Lagrangian is invariant under 0 —► — 0 and hence possesses a Z2 symmetry. 
For this reason, the potential has two minima: 0 = ±77, and the “vacuum 
manifold” has two-fold degeneracy. 

Consider the possibility that 0 = +7/ at z = +00 and 0 = —77 at 
z = —00. In this case, the continuous function 0 (z) has to go from —r] to 
+?7 as z is taken from —00 to +00 and so must necessarily pass through 
0 = 0 . But then there is energy in this field configuration since the potential 
is non-zero when 0 = 0 . Also, this configuration cannot relax to either of 
the two vacuum configurations, say 0(0) = +77, since that involves changing 
the field over an infinite volume from —77 to +77, which would cost an infinite 
amount of energy. 

Another way to see this is to notice the presence of a conserved current: 

77 = e^dycj) 

where 77, v = 0,1 and e * tv is the antisymmetric symbol in 2 dimensions. 
Clearly j 1 ' is conserved and so we have a conserved charge in the model: 

r/Q = J dzj° = 0(+oo) - 0(—00) . 

For the vacuum Q = 0 and for the configuration described above Q = 1 . 
So the configuration cannot relax into the vacuum - it is in a different 
topological sector. 

To get the field configuration with the boundary conditions 0 (±oo) = 
±77, one would have to solve the field equation resulting from the Lagrangian 
( 1 ). This would be a second order differential equation. Instead, one can 
use the clever method first derived by Bogomolnyi [ 2 ] and obtain a first 
order differential equation. The method uses the energy functional: 

E = J dz ^(<9 t 0) 2 + ^(9 2 0) 2 + K(0) 
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\{dtct>?+ \{d z <f> 
\{d t ct>) 2 + \{d z <t> 


V 2 V (<t>) f + V 2 v ^ d *4 

—— „i r<t>(+°°) 

2 V{4>) 


+ / d(j)'j2V {(/)') 

oo) V 


Then, for fixed values of 4> at Too, the energy is minimized if 


d t (t) = 0 


and 

dA - ^2V{<j>) = 0 . 

Furthermore, the minimum value of the energy is: 



In our case, 

y/vtf) = ^(7 ~ <t> 2 ) 

which can be inserted in the above equations to get the “kink” solution: 


0 = i] tanh 



for which the energy per unit area is: 


@kink 



1 m 3 

3 71 


( 2 ) 


where m = \/2Xr) is the mass of excitations (particles) of <f> in the vacuum 
of the model. Note that the energy density is localized in the region where 
(/> is not in the vacuum, i.e. in a region of thickness ~ m -1 around z = 0. 

We can extend the model in eq. (1) to 3+1 dimensions and consider the 
case when <f> only depends on z but not on x and y. We can still obtain 
the kink solution for every value of x and y and so the kink solution will 
describe a “domain wall” in the xy— plane. 

At the center of the kink, cf> = 0, and hence the Z 2 symmetry is restored 
in the core of the kink. In this sense, the kink is a “relic” of the symmetric 
phase of the system. If kinks were present in the universe today, their 
interiors would give us a glimpse of what the universe was like prior to the 
phase transition. 
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3. SU( 5) x Z 2 walls 

An example that is more relevant to cosmology is motivated by Grand 
Unification. Here we will consider the SU( 5) model: 

L = Tr {D^f - ^Tr(X^X^) - U($) (3) 

where, in terms of components, = 4> a T a is an SU(5) adjoint, the gauge 
field strengths are X^ = X^ u T a and the SU(5) generators T a are normal¬ 
ized such that Tr(T a T b ) = 5 ab / 2. The definition of the covariant derivative 
is: 

D^ = d^-ie[X il ^} (4) 

and the potential is the most general quartic in $: 

V($) = -m 2 TY($ 2 ) + h[ Tr ($ 2 )] 2 + ATr(4> 4 ) + 7 Tr($ 3 ) - V 0 , (5) 

where, Vo is a constant that we will choose so as to set the minimum value 
of the potential to zero. 

The SU( 5) symmetry is broken to 

H = [SU(3) x SU( 2) x U(l)]/Z 6 ( 6 ) 

if the Higgs acquires a VEV equal to 

$0 = ^=diag(2,2,2,-3,-3) (7) 

where 

= rn y = h+ —A . (8) 

For the potential to have its global minimum at $ = <J>o, the parameters 
are constrained to satisfy: 


A > 0 , A' > 0 . (9) 

For the global minimum to have U(4 > o) = 0, in eq. (5) we set 

V 0 = -jif. (10) 

The model in eq. (3) with the potential in eq. (5) does not have any 
topological domain walls because there are no broken discrete symmetries. 
In particular, the Z 2 symmetry under — 4> is absent due to the cubic 
term. However if 7 is small, there are walls connecting the two vacua related 
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by $ —► — $ that are almost topological. In our analysis we will set 7 = 0, 
in which case the symmetry of the model is SU( 5) x Z 2 . An expectation of 
breaks the Z 2 symmetry and leads to topological domain walls. 

The Lagrangian in eq. (3) provides us with the (second order) equations 
of motion for all the fields. However, it does not tell us the boundary con¬ 
ditions on the fields. The boundary conditions will depend on the class of 
solutions that we are interested in. This is the first issue that we need to 
settle. 

We would like a solution that has the correct topology. Topology tells us 
that the vacuum manifold is disconnected. One of the disconnected regions 
is described by —U^$qU and the other by -\-U^qU, where U € 577(5). 
To obtain solutions we need = $(— 00 ) to be in the first sector and 
4> + = 4>(+oo) to be in the second disconnected sector. Since we can globally 
rotate the fields by any SU( 5) transformation, we are free to set = —$o- 
But we still have to choose = +W$qU for some U £ 517(5). 

The freedom to choose is greatly reduced by the following result: a 
static solution to the field equations is only possible if 4> + ] = 0. I will 
not give the proof here but you can find it in the Appendix of Ref. [4]. 

The condition [$_,<I> + ] = 0 still permits a large variety of boundary 
conditions 4> + . But the different choices for <h + can all be rotated into a 
diagonal choice by a transformation in the unbroken global symmetry group 
which leaves <f>_ invariant. Hence we need only consider a discrete set of 
choices for 4> + , and these are: 


4>V 0) = e^diag(2,2,2,-3,-3) 




e^| di ag( 2 , 2 ,-3,2,-3) 

e^=diag(2, -3, -3,2, 2) 


( 11 ) 

( 12 ) 

(13) 


The superscript on 4> + corresponds to the number of entries of that 
have been permuted. The parameter e has the value —1 for the topological 
domain walls. Later we will find that non-trivial solutions also exist for 
e = +1. These correspond to non-topological domain walls. 

It is tempting to try to obtain the domain wall solutions using Bogo- 
molnyi’s method. This has not been done and remains an interesting open 
problem. However, we can write down the solution for the zero permutation 
case since it is exactly like the Z 2 case: 


= tanhpr)$ 0 


(14) 


where X = mxjyj2. 
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The other cases are a little more involved but still quite simple. The 
solution is of the form: 


= F + (x)M + + F_(x)M_ + g(x)M. , (15) 


where 


M + 


4> + + <f>_ 
2 


M_ 



(16) 


g(±oo) = 0 and M is yet to be found. 

It is now convenient to move from the SU (5) to the SU ( N ) case. Let 
N = 2n + 1, so n = 2 when N = 5. Also let us label the number of 
permutations in the boundary conditions by the integer q. For N = 5, the 
three possible boundary conditions in eq. (11)-(13) correspond to q = 0,1, 2. 
The advantage of dealing with the general case is that we can write one 
equation in terms of n and q rather than 3 equations, one per case. 

The formulae for M± can now be explicitly written (we set e = — 1 for 
now): 

M+ = //A f ^/^^|^J^^ydiag(O n _|_l_g, l g , —lg, On-g) , (17) 

M— = _Dydiag( — 2nl n+ i_ g , 1 9 , l q , 2 (n + l)l n _ g ) . (18) 


We have used 0^ and 1 j. to denote the k —dimensional zero and unit matrices 
respectively. Note that the matrices M± are orthogonal: 


Tr(M + M_) = 0 , (19) 

but are not normalized to 1/2. The boundary conditions for F± are: 

F-(— oo) = -1 , F_(+oo) = +1 , 

F + (—oo) = +1 , .F_|_(+oo) = +1 . (20) 


The advantage of this form of the ansatz is that, for particular values of the 
parameters of a quartic potential in the q = n topological (e = — 1) case, 
one finds the explicit and simple solution F + (x) = 1, Tb (x) = tanh(crx) 
and g(x) = 0, where a = mx/y/2 [5, 4], Also, for q = 0, e = — 1, the solution 
is the embedded Z 2 kink i.e. F + (x) = g{x) = 0, F_(x) = tanh(crx). 

Now we would like to find the unknown matrix M in the ansatz (15). 
By inserting the ansatz in the equations of motion [6], it turns out that M 
is uniquely determined and this is: 


M = /i diag ( q(n - q)l n+ i- q ,-(n - q)(n + 1 - q)l q , 

- {n-q){n + l-q)l q ,q{n + l-q)l n - q ) (21) 
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with /i being a normalization factor in which we also include the energy 
scale r] for convenience: 




= r)[2q(n — q)(n + 1 — g){ 2 n(n + 1 — q) — 7 }] 1 / 2 . 


( 22 ) 


Note that the matrix M is not normalizable if q = 0 or if q = n. For these 
values of q, we can set g{x ) = 0 . 

The functions F±(x ) and g(x) can be found by solving their equations of 
motion derived from the Lagrangian together with the specified boundary 
conditions. There is no guarantee that a solution will exist and so we find 
the solutions explicitly for N = 5 with a quartic potential. The solution for 
q = 0 is simply that in eq. (14). The <7 = 1 profile functions are evaluated 
numerically and shown in Figure 1. The profiles for q = 2 have also been 
evaluated numerically in Ref. [5]. However, an interesting situation arises 
for q = 2 when 

h 3 3 (23) 


A N(N — 1) 
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For this special value of coupling constants, the profile functions can be 
written as [4]: 

F + = 1 , F_ = tanh(AC) , g = 0 (24) 

It is not understood what is special about h/X = —3/20 that it allows an 
analytic solution. Could it be that there is a symmetry of the system at 
this point in parameter space that simplifies matters? Could a Bogomolnyi 
type analysis be done at this point? 

The stability of these solutions has been investigated in Ref. [5, 4, 6 ] 
and the q = 2 solution turns out to be the minimum energy domain wall 
solution. For the particular choice of model parameters given in eq. (23) 
(valid for N > 3) [4], the energy of the q = 2 kink is: 


o = 


2\/2 


m 


3 A \N — 3 


N - 1 


(25) 


An interesting point to note is that the ansatz is valid even if T-t are not 
in distinct topological sectors i.e. even if e = +1. These imply the existence 
of non-topological kink solutions in the model. If we include a subscript 
NT to denote “non-topological” and T to denote “topological”, we have 

&NTk = F + (x)~Mnt+ + F 1 _(x)MtvT- + g(x)M.NT ■ (26) 

Since "Fatth- = — 4>t+> we find 


Mjvt+ = M T — , ]VIjv2 ’_ = Mr + , Mjvt = . 


(27) 
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Figure 1. The profile functions F+ (nearly 1 throughout), F- (shaped like a tanh 
function), and g (nearly zero) for the q = 1 topological kink with parameters h = —3/70, 
A = 1 and rj = 1. 


Hence 

3? NTk = F_(x)Mt+ + F+(x) Mr_ + g(x)~M.T • (28) 

So to get F_ (F + ) for the non-topological kink we have to solve the topo¬ 
logical F + (F_) equation of motion with the boundary conditions for F- 
(F + ). To obtain g for the non-topological kink, we need to interchange F + 
and F- in the topological equation of motion. The boundary conditions 
for g are unchanged. Generally the non-topological solutions, if they exist, 
will be unstable. However, the possibility that some of them may be locally 
stable for certain potentials cannot be excluded. 

3.1. DOMAIN WALL LATTICE 

The forces between domain walls of different kinds has been studied by 
Pogosian [7]. What is most interesting is that the q = 2 wall and anti-wall 
can repel. Here one needs to be careful about the meaning of an “anti-wall”. 
An anti-wall should have a topological charge that is opposite to that of 
a wall. That is, a wall and its anti-wall together should be in the trivial 
topological sector. But this condition still leaves open several different kinds 
of anti-walls for a given wall. To be specific consider a domain wall across 
which the Higgs field changes from (2, 2, 2, —3, —3) to —(2, —3, —3, 2,2). (I 
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am suppressing the normalization factor for convenience of writing.) There 
can be two corresponding anti-walls. In the first type - called Type I - the 
Higgs field can go from — (2, —3, —3,2, 2) to +(2, 2,2, —3, —3), thus revert¬ 
ing the change across the first domain wall. In the second type (Type II), 
the Higgs held can go from —(2, —3, —3, 2, 2) to +(—3, 2, 2, —3, 2). Pogosian 
found that the force between a wall and its Type I anti-wall is attractive, 
but the force between a wall and its Type II anti-wall is repulsive. This 
is understood by noting that the force between walls is proportional to 
T r (Qi£? 2 ) where Q\ and Q 2 are the charge matrices of the two walls. If 
the Higgs held to the left side of the walls is 4+, between the two walls is 
$ 0 ) and is on the right-hand side of the walls, then Q\ oc d>o — 4+ and 
Q 2 oc — $ 0 - Then the stable (q = 2) walls can have charge matrices 
proportional to (-4,1,1,1,1), (1,-4,1,1,1), (1,1,-4,1,1), (1,1,1,-4,1) 
and (1,1,1,1, —4). (Stable anti-walls have the same charges but with a 
minus sign.) Then it is easy to see that one can take a wall with one of 
the 5 charges listed above and it will repel an anti-wall which has the —4 
occurring in some other entry. Hence there are combinations of walls and 
anti-walls for which the interaction is repulsive. 

We know that topology forces a wall to be followed by an anti-wall. 
Then we can set up a sequence of walls and anti-walls in the following way: 

...Q (1) Q (5) Q (3) Q (1) Q (5) Q (3) ... (29) 

where Qi and Qi refer to a wall and an anti-wall of type i respectively. 
Alternately, this sequence of walls would be achieved with the following 
sequence of Higgs field expectation values: 


-(2, 2, 2,-3,-3) 


+(2,-3,-3, 2, 2) 
-(-3, 2, 2,-3, 2) 
+(2,-3, 2, 2,-3) 
-(2, 2,-3,-3, 2) 
+(-3,-3, 2, 2, 2) 
-(2, 2, 2,-3,-3) 


(30) 


The forces between walls fall off exponentially fast and hence the dominant 
forces will be between nearest neighbors. 

Note that the sequence of walls is periodic with a period of 6 walls, and 
these 6 walls have a net topological charge that vanishes. Hence we can put 
the sequence in a periodic box i.e. compact space. This gives us a finite 
lattice of domain walls. 

The sequence described above is one whose period is the minimum pos¬ 
sible (namely, 6). It is easy to construct other sequences with greater peri- 
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odicity. For example: 

...Q (1) Q (5) Q (3) Q (4) Q (2) Q (1) Q (5) Q (3) Q (4) Q (2) ... (31) 

is a repeating sequence of 10 domain walls. 

In ongoing work, we have checked that the lattice is a solution of the 
equations of motion. However, it is not stable. The instability is towards 
rotations of the walls in internal space. For example, the Q ^ antiwall could 
rotate to become either a Q ^ or a Q 1 ' ^ antiwall. Once this rotation takes 
place, this antiwall will be attracted by the Q ^ or the neighbor, and 
annihilation can occur. An analogy of this process is due to Pogosian - 
if we place several bar magnets along a line such that the North poles of 
neighboring magnets face each other, then each bar magnet will be unstable 
to rotation. Once the bar magnets rotate and North poles face neighboring 
South poles, the magnets will attract and the chain of magnets is unstable 
to collapse. In this analogy, the rotation is in physical space; in the case of 
domain walls, the rotation is in internal space. 

The understanding of the domain wall lattice instability, provides us 
with another avenue of investigation. If the bar magnets were confined to 
one dimensions (for example, they could be placed in a tube, or attached to 
a wire), then the instability would cease to exist since rotations would not 
be possible. The corresponding situation in the case of domain walls is if the 
unbroken symmetry is too restrictive to allow rotation. This can be achieved 
by breaking the SU(5) x Z -2 to U{ l) 4 or something smaller. A simpler 
alternative is to start out with a model in which there are no rotational 
degrees of freedom present. This scheme is discussed in [8]. One starts 
with a model that only has a discrete symmetry group that corresponds to 
permutations of the diagonal entries of SU(5 ) and the sign flip given by 
Z- 2 - It can be shown that the domain wall lattice is stable in this model. 

3.2. FORMATION OF DOMAIN WALLS 

The properties of the network of domain walls at formation has been deter¬ 
mined by numerical simulations. The idea behind the simulations is that the 
vacuum in any correlated region of space is determined at random. Then, 
if there are only two degenerate vacuua (call them + and —), there will be 
spatial regions that will be in the + phase and others in the — phase. The 
boundaries between these regions of different phases is the location of the 
domain walls. This is nothing other than the “Kibble mechanism”. 

By performing numerical simulations, the statistics shown in Table I was 
obtained [9]. The data shows that there is essentially one giant connected + 
cluster. By symmetry there will be one connected — cluster. In the infinite 
volume limit, these clusters will also be infinite and their surface areas will 
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Cluster size 

1 

2 

3 

4 

6 

10 

31082 

Number 

462 

84 

14 

13 

1 

1 

1 


also be infinite. Therefore the topological domain wall formed at the phase 
transition will be infinite. 

What does the Kibble mechanism predict for SU (5) x Z 2 domain walls? 
Once again we have to throw down values of the Higgs field on a lattice and 
then examine the walls that would form at the interface. We have found 
that there are 3 kinds of wall solutions and so each one will be formed 
with some probability. In the Kibble mechanism approach, the probability 
that a certain wall will form is directly related to the number of boundary 
values that result in the formation of a defect. So we need to evaluate all 
the boundary conditions that will lead to domain walls with a certain value 
of q. In other words, we want to determine the “space of kinks” for a fixed 
value of q [4, 6]. 

Let the Higgs field to the left and right of a domain wall be and 4># 
respectively. Without loss of generalization, we can take = To- To start, 
let us find all values of Tr that will give a q = 0 kink. But there is only 
one such value, namely = — ‘ho- So the space of the q = 0 kink is just 
one point. Next, let us determine all 4*^ that will give a q = 1 kink. Here 
there are several possibilities since different entries of 4>o can be permuted 
to give different In fact, one can take any choice of for example 
along —(2, 2, —3, 2, —3) and act by gauge transformations belonging to the 
unbroken symmetry group H (see eq. (6)), and this will lead to another 
choice for 4>^. Of course, some of these rotations will leave unchanged 
and these should not be counted. So the space of q = 1 kinks is given by 
the coset space: 

Hi = H/K\ (32) 

where 

K\ = [SU{ 2) x U(lf}/Z 2 . (33) 

A similar argument shows that the space of q = 2 kinks is given by: 

H 2 = H/K 2 , (34) 

where 

K 2 = [SU{2f x U(lf]/Zl . (35) 

Now that we have the space of kinks for each q, we note that the space 
of the q = 0 kink is zero dimensional, the space of the q = 1 kink is 6 
dimensional, and the space of the q = 2 kink is 4 dimensional. Hence, in 
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the Kibble mechanism approach, the probability of a kink being of the q = 0 
or q = 2 variety is zero, and the probability of the kink being of the q = 1 
variety is 1. 

A subtlety that has not been discussed above is that there is also the pos¬ 
sibility that if we lay down Higgs fields randomly, we may get [<J>£, <h/j] ^ 0. 
In this case, as described earlier, there will be no static solution to the 
equations. Then the field configuration will evolve towards a static config¬ 
uration. Our discussion above assumes that such a configuration has been 
reached, and neighboring domains always have values of d> that commute. 
This is not completely satisfactory since there will be time scales that are 
associated with the relaxation and these must be compared to other time 
scales characterizing the phase transition. 

To summarize, the Kibble mechanism predicts that only q = 1 domain 
walls will be formed at the SU (5) x Z 2 phase transition. However, we know 
that the stable variety of walls have q = 2, and the q = 1 walls will decay 
into them. Precisely how the q = 2 walls convert into q = 1 walls during a 
phase transformation has not been studied. Nor is it known how long this 
relaxation will take. The answers to these questions is of some importance 
as we will briefly see in the next section. 

3.3. IMPORTANCE IN COSMOLOGY 

The SU (5) x Z 2 symmetry breaking also leads to topological magnetic 
monopoles which are cosmologically unacceptable. However these magnetic 
monopoles interact with the domain walls that we have discussed in the 
earlier sections. If the domain wall is of the q = 0 variety, the monopole is 
certainly destroyed and its magnetic charge spreads out on the domain wall 
[3, 5]. If the domain wall is of the q = 1 variety, monopoles in a large number 
of orientations in the group space get destroyed whereas others survive. If 
the domain wall is of the q = 2 variety, a smaller number of monopoles 
get destroyed [6]. So the cosmological predictions depend on the kinds of 
walls that are formed and their interactions with magnetic monopoles. This 
is why it is important to consider the types of domain walls that will be 
formed and the eventual relaxation of the system, especially in regard to 
the number density of magnetic monopoles. 
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